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We present a study of transport in graphene devices on polar insulating substrates by solving 
the Bolzmann transport equation in the presence of graphene phonon, surface polar phonon, and 
Coulomb charged impurity scattering. The value of the saturated velocity shows very weak depen- 
dence on the carrier density, the nature of the insulating substrate, and the low-field mobility, varied 
by the charged impurity concentration. The saturated velocity of 4 - 8 x 10 7 cm/s calculated at 
room temperature is significantly larger than reported experimental values. The discrepancy is due 
to the self-heating effect which lowers substantially the value of the saturated velocity. We predict 
that by reducing the insulator oxide thickness, which limits the thermal conductance, the saturated 
currents can be significantly enhanced. We also calculate the surface polar phonon contribution to 
the low-field mobility as a function of carrier density, temperature, and distance from the substrate. 

PACS numbers: 72.80. Vp, 72.10.Di, 73.50.Fq 



I. INTRODUCTION 

The excellent transport^— and optical properties of 
graphene^ have attracted strong interest in possible ap- 
plications of this material in nanoscale electronics and 
optoelectronics^—. The electrostatic modulation of the 
graphene channel through gates yields very promising 
two-dimensional field-effect devices for analog and radio- 
frequency application ;) 9 ' 10 . Such devices should ideally 
be operated in the saturation limifeii. Indeed, it has 
been shown that the current saturates as the source-drain 
field is increased to a few Volts per micronii - — . While 
clastic scattering determines the low-field mobility, the 
current saturation has been attributed due to the inelas- 
tic scattering by either surface polar phonons (SPP) of 
the polar substrate s 11 ' 12 or the intrinsic graphene optical 
phonons^. In addition to the uncertainty on the nature 
of the inelastic scattering mechanism, significant heating 
of the graphene devices operated under high bias con- 
ditions is expected. This has recently been measured by 
Raman spectroscop y 12 ' 16 . However, little is known about 
the role of self-heating and elastic scattering on the cur- 
rent saturation. 

High bias measurements in grapheneii~— were ana- 
lyzed by analogy to ID carbon nanotubes, where the 
magnitude of the saturated current is determined by 
the optical phonon energy responsible for the saturation. 
However, a priori an extension of the simple analytical 
model for the saturated velocity in 2D graphene, as an 
inverse of the square root of the carrier density, may not 
be applicable. In RefJ^, using hydrodynamic transport 
theory including graphene optical phonon scattering, the 
saturation velocity was found to be weakly carrier den- 
sity dependent, whereas in RefJ^, using a Monte Carlo 
solution of the Boltzmann transport equation (BTE) , the 
saturated velocity was found to follow an inverse square 
root dependence on the carrier density. 

In this work, we explore the effects of Coulomb impu- 
rity, graphene phonon, and SPP scattering on different 
polar substrates on the current saturation in the diffu- 



sive transport regime. We find that the self-heating of 
graphene on Si02 limits significantly the value of the 
saturated current. The electronic structure of graphene 
is described by a 7r-orbital tight-binding model with a 
hoping parameter to = 3.1 cV, which gives a Fermi ve- 
locity v F = (V3/2)t a/h fn 10 6 m/s, where a = 0.246 
nm is the graphene lattice constant. For the electron- 
phonon scattering we use the Su-Schrieffer-Heeger (SSH) 
modeU£ to express the modulation of the 7r-orbital over- 
lap t = t(j—g8Rcc with C-C distance Rcc- The electron- 
optical phonon scattering has been calculated by first 
principle s 18 ' 19 . It was found that SSH reproduces fairly 
well the LDA results for the electron-optical phonon scat- 
tering if g sa 4.5 eV/Ai£. On the other hand, electron- 
electron correlations, taken into account using the GW 
approximation 2 ^, give electron-K-point optical phonon 
coupling corresponding to g ~ 6.5 cV/A and coupling to 
T-point phonons to give g w 5.3 eV/A. Therefore, it is ex- 
pected that electron-phonon coupling in graphene can be 
renormalizcd depending on the environment which will 
screen the electron-electron interactions. In this work we 
use an average value of g = 5.3 eV/A as in Refi 2 ^. The 
parameters used for SPP scattering on Si02, Hf02, SiC, 
and BN polar substrates are given in Tabic Q] 



II. LOW-FIELD MOBILITY 

The low-field mobility in pristine graphene, in the ab- 
sence of charged impurities and defects, is determined by 
scattering from the graphene phonons and it is shown 
in Fig. Q] There are two acoustic phonon branches, 
transverse (TA) and longitudinal (LA), with an ap- 
preciable electron-phonon coupling. Within the SSH 
model, the TA and LA modes have different angle depen- 
dencies for the electron-phonon couplings \Mk k+q\ 2 = 
D 2 ac q 2 h/{8NM C io q )(l ± cos{3{e k + 6 k+q )))2Z, where D ac 
is a deformation potential, 9 k is a directional angle of 
wavevector A»21, u) q is a phonon frequency, Mc is a mass 
of carbon atom, and A is a number of k-points. There- 
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fore, the acoustic phonon coupling can be approximated 
by \M k , k+q \ 2 = Dl c qh/(4NM c v ph ), where v ph is a char- 
acteristic (LA/TA) sound velocity. The value of deforma- 
tion potential is given by2£: D ac = 3gan/ (4-^/3), where 
the reduction factor k = v\ A j (v\ A — v\ A ) was introduced 
in Refi 28 ' 29 . For the valence phonon models used here 
we obtain D ac ss 3.7 eV, consistent with the numerical 
calculations. 

The low-field mobility [i can be found using Boltzmann 
theory a — enfi = e 2 v F D n T /2, where a is the conduc- 
tivity, n is the carrier density, D n — 2Ep / (nh 2 vp) is the 
density of states, Ep ss hvpy/irn is the Fermi energy, and 
r is the scattering time. The latter can be found as: 

- = IT y2\M k ,k +q \ 2 (l-cos(8 k -0 k+q )) 

i 

(N q 6(E k - E k+q + hw q ) + 

(N q + l)S(E k - E k+q - hw q )) (1) 

where N q is the Bosc-Einstein phonon occupation num- 
ber. The summation in Eq. ([1} is replaced by the integral 
(1/N) J2 q = A/ (47r 2 ) J qdqdO (sum over one spin and one 

valley), where A = \/3a 2 /2 is the area of the two atom 
unit cell. The low-field mobility and the scattering rate 
in the high temperature T limit are given by: 

ep m ftv 2 F v 2 h l 
1 1 E k D 2 ac 

— = i3 72 —kbT (2) 

where p m = 2M C /A m 7.66 x 10~ n kg/cm" 2 is the 
graphcnc mass density. Note a numerical difference of a 
factor of 4 in Eq. © for r used in the literatur e 26 ' 31 . This 
discrepancy leads, in part, to the large range of deforma- 
tion potential values from D ac =10 eV to 30 cV quoted in 
the literature 3 -^— . In addition, uncertainty in the sound 
velocity can also contribute to the spread of the values of 
the deformation potentials. For example, the deforma- 
tion potential extracted from the resistivity temperature 



dependence in graphcnc on SiC>2 at n — 10 12 cm 2 in 
RefM would give D ac = 7.1 cV, if Eq. © is used with 
Vph = 17.3 km/s22. The measurements for the electron 
branch in suspended graphene at n = 2 x 10 11 cm~ 2 in 
Refi 3 ^ would correspondingly give D ac = 12 eV. When 
quantifying the acoustic scattering in different studies it 
would, therefore, help to report not only the value of D ac , 
but also that of the low-field conductivity extrapolated 
to room temperature. In our model, it corresponds to 
cr ac (T = 300AT) fn 0.1 S. Note that in low mobility sam- 
ples, due to charged sites in the substrate, an additional 
temperature dependence from Coulomb scattering may 
ariseai. 
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TABLE I: Parameters for the SPP scattering for graphene 
on Si02, Hf02, SiC, and hexagonal BN substrates. The 
surface optical phonon (SO) energies are obtained from 
the bulk longitudinal optical (LO) phonons as huso — 





SiO^ 


HfO^ 


SiG^ 


h-BN2£ 


eo 


3.9 


22.0 


9.7 


5.09 


Ei 


3.36 


6.58 




4.575 


Eoo 


2.40 


5.03 


6.5 


4.10 


hwsoi in meV 


58.9 


21.6 


116.0 


101.7 


hwso2 in meV 


156.4 


54.2 




195.7 


Fx in meV 


0.237 


0.304 


0.735 


0.258 


Ff in meV 


1.612 


0.293 




0.520 



FIG. 1: (Color online) Temperature dependence of the low- 
field mobility in graphene on different substrates from top 
to bottom: green circles - intrinsic, black open circles - on 
BN, cyan squares - on SiC, red triangles - on Si02, blue open 
diamonds - on Hf02 at carrier densities (a) n = 10 12 cm~ 2 
and (b) n = 5x 10 12 cm' 2 . The solid curves are results using 
Eq. ([6]) with characteristic sound velocity v p h = 17.3 km/s, 
energy of the optical phonon in graphene hu) ov = 0.19 eV. 



The two optical phonons at the T point have couplings 
\ M ki+ q \ 2 = D 2 h/(2NM c uj r )(l ± ss'cos(9 k - 6 k+q )) for 
LO (- sign) and TO (+ sign) modes respectivel y 18 ' 38 , 
where D T = 3g/V2 w 11.2 cV/A 3 ^, s = 1 for elec- 
trons and s = — 1 for holes. The K-point TO phonon 
mode has an electron-phonon coupling twice as larg o 18 ' 21 , 
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FIG. 2: (Color online) Carrier density of the low-field mo- 
bility in graphene on HfC>2 at T=300 K for distances from 
the substrate (from top to bottom) zo = 12.5 A (cyan open 
squares), 10.0 A (red open circles), 7.5 A (magenta triangles), 
5.0 A (blue squares), 3.5 A (green circles). The solid curves 
are results using Eq. ©. 
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with the angle dependence given by 
D*h/(NMcu K )(l - ss'cos(6 k - 9 k+q ))^. The effect of 
the optical phonons (both at V and K) on the low-field 
mobility can be calculated according toM: 



ep m v F uj p 



1 



2irDl p 



nN n 



(3) 



where D op = 2Dr = 22.4 eV/A is an effective electron- 
optical phonon coupling. The angle integration in Eq. ([T]) 
for A"-point phonons gives a factor of 3/2 larger contri- 
bution than that for T-point phonons scattering. 

The SPP phonons on polar substrates produce an elec- 
tric field that couples to the electrons on the nearby 
graphene. While the SPP scattering is of lesser impor- 
tance in conventional MetalOxidcScmiconductor Field- 
Effect Transistors (MOSFETs} 2 ^,, it is much more promi- 
nent in graphene and carbon nanotubes due to the much 
smaller vertical dimension of the devices given by the 
van der Waals distance. It has been invoked to explain 
transport under both low and high bias conditions in 
graphen o 11 ' 35 ' 41 ' 42 and carbon nanotubes 4 ^ on polar sub- 
strates. In graphene the SPP coupling is given b y 41 ' 42 : 



l + ss'cos(6 k -6 k+q )4Tr 2 e 2 F2 
2 NAq 



-2qz 



(4) 



where z ss 3.5 A is the van der Waals distance between 
the polar substrate and the graphene flake. The mag- 
nitude of the polarization field is given by the Frohlich 



coupling^: F 2 = ( eoo + Eenv ~ ^R^) . where 

fouso,v is a surface phonon energy and Sq and Soa are the 
low- and high-frequency dielectric constants of the polar 
substrate, see Table |U The screening of the Coulomb 
interaction by the environment above the polar dielec- 
tric is taken into account by e cnv . Since the screening of 
the electric field perpendicular to the graphene plane is 
weak 4 ^, we take e cnv = 1. When there are several SPP 
modes with an appreciable coupling, then the low- and 
high- frequency e are understood as intermediate dielec- 
tric functions at uji <C ^so,v for eo and at w% 3> ojso,v 
for eoJZ. 
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FIG. 3: (Color online) Current density electric field depen- 
dence in graphene on different polar substrates: green circles 
- intrinsic, black open circles - on BN, cyan squares - on SiC, 
red triangles - on Si02, blue diamonds - on Hf02 at carrier 
densities (a) n = 10 12 cm -2 , (b) n = 5 x 10 12 cm -2 , and (c) 



n = 10 cm and a lattice temperature T am b 
solid curves are fits to Eq. (0. 



300 K. The 



As shown in Fig. [T] the SPP contribution to the low- 
field mobility can be well approximated by: 



P'Spp.u 



f3 Hvf e,vp exp (fco z o) 
hu v e 2 F 2 N spp 



(5) 



which is a non-monotonic function of carrier density 
n. Here fco ~ \J {2loso.u/ v f) 2 + ocn, where parameters 
a w 10.5 and (3 rj 0.153 x 10~ 4 eV are determined to 
give an overall agreement of the mobility dependence 
on the carrier density n and distance z , as shown in 
Fig. [2J SiC grown graphene has an intermediate dead 
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layer—, which increases the effective distance zq and also 
changes the environmental screening e env . Both effects, 
have to be taken into account when direct comparison 
is made with the experiment. Note that the calculated 
temperature dependencies shown in Fig. [1] deviate from 
the scattering rate temperature dependence given by the 
SPP phonon occupation number N spPlV . An additional 
temperature dependence arises from the thermal averag- 
ing of the scattering rate with the carrier distribution 
function. We have chosen a set of parameters a and (3 to 
agree best with the mobility at room temperature. 

The calculated low-field mobility from the BTE solu- 
tion in Fig. Q] can be well described using Matthicsscn's 
rule: 

M _1 = Mac 1 + Hop + ^VVV ( 6 ) 
a u 

where mobility contributions due to the acoustic, optical, 
and SPP phonons are given by Eq. @, and (|5|), 
respectively. 
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III. CURRENT SATURATION 

At high bias, the transport in graphene is usually de- 
scribed by the saturated current mode l 11 ' 21 ' 47 : 



where r\ = n e + rih is the total carrier density, v sa t is 
a saturation velocity, and F is an electric field. Note 
that in the model ^ e = Hh- Within simple model, in 
the full saturation regime only carriers around the Fermi 
energy Ep in the energy window Ep ± Ml/ 2 contribute 
to the current. Here Vl is the characteristic frequency of 
the phonon responsible for the current saturation. The 
saturated velocity can then be calculated and for Ep > 
Ml/2: 

2hv F 2 n 

v sa t ~ — = -= (8) 

7r Ep ir y/irn 

Note again a difference in the numerical prefactor used 
in the literature^— 1 versus 2/ir which affects the nu- 
merical value of the characteristic phonon energy Ml ex- 
tracted from the experimcnt a 11 ' 12 . 
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FIG. 4: (Color online) (a) Current saturation in graphene on 
S1O2 at different temperatures T in K and Coulomb impu- 
rity concentration n; in 10 12 cm -2 : (T,rii)= (300, 0.0) - red 
circles, (300, 1.35) - red squares, (600, 0.0) - blue triangles, 
(600, 2.8) - blue diamonds. The solid curves are the best 
fits to Eq. (0) and the corresponding saturated velocities are 
v a at = 7.26, 7.07, 7.15, 6.91 xlO 7 cm/s, which are the same 
values within 5%, despite the fact that the low- field mobilities 
vary by a factor of 15. (b) calculated total carrier density as 
a function of the electric field for the same cases as in (a). 
The charge carrier density was fixed at n = n e — «/,= 2 x 10 12 
cm -2 . 



FIG. 5: (Color online) (a) Reference velocity v re f ~ v sa t/2 
and (b) reference field F re f ~ v sa t / /i dependencies on the low- 
field mobility. The results are obtained from the calculations 
at room temperature for graphene on Si02 for different impu- 
rity concentrations, which determine the low-field mobility^. 
Calculations are for fixed charge carrier densities n (in 10 12 
cm -2 ) of 1.0 - red circles, 2.0 - blue squares, 3.0 - green di- 
amonds, 5.0 - cyan triangles. The black dashed curve in (b) 
shows F re f = Vsat/fJ. dependence with v sa t = 8 x 10 7 cm/s. 

The current densities for graphene on different sub- 
strates as a function of electric field are shown in Fig. [3] 
When phonons are kept in thermal equilibrium at T amo = 
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300 K, the current docs not show full saturation for the 
experimentally relevant source-drain fields up to 2 V//im. 
In Fig. |3k, the current shows negative differential con- 
ductance for scattering by intrinsic graphene phonons. 
This behavior is due to the deviation of the bandstruc- 
ture from the linear band dispersion, similar to the effect 
of the non-parabolicity in carbon nanotubes^i. While 
the current at high bias in Fig. [3] is found to be larger 
for graphene on polar insulators with larger SPP phonon 
energy, the saturated velocity as obtained from the fit to 
Eq. (|7J| does not obey Eq. ©. 
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FIG. 6: (Color online) Carrier density dependence of the sat- 
urated velocity in graphene on different substrates at room 
temperature: green circles - intrinsic, black open circles - on 
BN, cyan squares - on SiC, red triangles - on Si02, blue di- 
amonds - on Hf02- The dashed curves are predictions using 
Eq. (8J with optical phonon energy hwsoi from Tabled] from 
top to bottom: green - intrinsic, cyan - on SiC, black - on BN, 
red - on Si02, blue - on Hf02- 

First we explore the dependence of the saturated ve- 
locity on the low-field mobility. The latter is strongly 
affected by the quality of the substrate, which in- 
troduces Coulomb impurity scatterin g 35 ' 37 ' 48 " — , reso- 
nant scatterers 5 ^— , electron-hole puddles 5 -^, and surface 
roughnes a 58 ' 59 . In Fig. 0^ we show current- voltage char- 
acteristics of graphene on Si02 with a variable Coulomb 
impurity concentration rii following Refi^i. Despite the 
fact that the low- field mobilities vary by a factor of 15, 
the saturated velocity, as obtained form the fit to Eq. (J7J , 
is essentially the same. The electronic temperature in- 
creases with the electric field, such that minority carriers 
contribution to the current becomes significant. This ef- 
fect is more prounced at elevated lattice temperatures as 
shown in Fig. 2Jd. 

The systematic behavior of the low-field mobility, 
which is modified by the impurity concentration, on the 
saturated velocity is shown in Fig. [5J Since the current 
does not fully saturate even at large fields, we have cho- 
sen to plot in Fig. [S] a reference velocity v re f calculated 



at the reference field F re f where the mobility is reduced 
by a factor of two. If Eq. is applicable, then the 
reference velocity is equal to half the saturated velocity 
v sa t and is reached at a reference field of F re f = v sa t/fJ-- 
We find that the reference velocity is about 3 — 5 x 10 7 
cm/s and has little dependence on either carrier density 
or low-field mobility in the parameter space considered 
here. The reference field is inversely proportional to the 
low-field mobility as shown in Fig. [5}), which justifies the 
functional form Eq. Q used to analyze the BTE results. 



IV. SELF-HEATING EFFECT ON CURRENT 
SATURATION 

The analysis of both experimental^— and 
simulation^ 5 - results had often relied on Eq. (|SJ). In 
Fig. [5] we show that Eq. ((5J) fails qualitatively to describe 
the results of the BTE simulations. In particular, 
saturated velocity in graphene on Si02 is predicted to 
be very similar to that in graphene on BN and SiC 
substrates, although SPP phonon energies in the latter 
are almost a factor of two larger. While Eq. © predicts 
a variation of v sa t by more than an order of magnitude 
for the range of n and Ml used in Fig. [SJ we find the 
values of v sa t from the fits to Eq. {7} to be within 
4 — 8 x 10 7 cm/s. The saturated velocity in intrinsic 
graphene shown by green curve in Fig. [6] agrees well 
with the Monte Carlo BTE solution in low density limit 
in Refi^S. At the same time, the reported experimental 
value o 11 ' 12 of v sa t at densities of 10 13 cm -2 are below 
10 7 cm/s, which is at least a factor of 4 smaller than 
that in Fig. [3 

We suggest that one of the factors for the apparent 
discrepancy is the self- heating effect. The temperature 
rise was shown^ to be proportional to the Joule losses 
T = T am b + jF/r, where T am b is the ambient tempera- 
ture and the thermal conductance r controls the heat dis- 
sipation. This simple picture applies only for graphene 
regions away from the contacts such that both the ther- 
mal contact resistance and the substrate thermal con- 
ductivity determine the value of r. Note that the upper 
bound for r, which corresponds to zero contact thermal 
resistance, is determined by the thermal conductivity and 
thickness of the insulating substrate. For example, ther- 
mal conductivity of Si02 is k w 1.4 W/(mK)£i and for 
the insulator height of h = 300 nm a maximum value 
of r = K/h « 0.47 kW/(K cm 2 ) is expected. However, 
this upper bound can be significantly reduced due to the 
thermal contact resistance as in Rcf. 12 . 

In the presence of the SPP scattering the role of the 
thermal contact resistance is minimized because electrons 
can give some of their energy directly to the substrate 
SPP phonon o 62 ' 63 . In our "self-heating model" we as- 
sume that SPP and graphene phonons are heated to the 
same temperature, which is proportional to the Joule 
losses found self-consistcntly, i.e. j = j(T). In sub- 
strates with much higher thermal conductivities, such as 
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FIG. 7: (Color online) Current density electric field depen- 
dence in graphene on different polar substrates: green circles 
- intrinsic, black open circles - on BN, cyan squares - on SiC, 
red triangles - on Si02, blue diamonds - on Hf02 at car- 
rier densities (a) n — 10 12 cm~ 2 , (b) n = 5 x 10 12 cm~ 2 , 
and (c) n = 10 13 cm -2 at self consistent lattice temperatures 
T = T amb + j(T)F/r for r = 0.47 kW/(K cm 2 ), correspond- 
ing to the ideal thermal contact resistance with 300 nm Si02- 
The vertical axis on the right shows corresponding values of 
the maximum temperatures at 2 V/^tm. 



the graphene phonon and SPP phonon energies. Indeed, 
the current density in the diffusive regime is j oc Ft. Ac- 
cording to Eq. flTJ), the self- heating reduces the scattering 
time as roc (1 + 2N op )~ 1 , such that j oc F/(l + 2N op ). In 
the high temperature limit N op oc ksT /fttv op and current 
saturates as j oc F / (fcBT am (,+ j F/r) , where we have used 
the self-heating temperature fc^T = kBT arn b+jF/r. The 
saturation current in this simple picture depends on r as 
j sa t oc r 5 , which is in qualitative agreement with a full 
BTE solution shown in Fig. \E]p. On the other hand, the 
large intrinsic graphene optical phonon leads to an ac- 
tivated temperature dependence of the scattering time, 
which produces the negative differential conductance in 

Fig.ra 
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SiC, BN and Hf02, the thermal contact resistance would 
dominate the value of r. In this case a full self-consistent 
solution including thermal transport in the substrate, 
which determines the SPP phonon temperature, would 
be needed. In Fig. [7J we find that the current densities 
drop by up to a factor of four at high biases as a result of 
self-heating for r = 0.47 kW/(Kcm 2 ). Moreover, the cur- 
rent now shows true saturation at experimental source- 
drain fields, which are much smaller than F Ie f calculated 
in Fig. [5Jd at room temperatures, and the high bias cur- 
rents are comparable to those reported in Refi 11 ' 12 . 

At high density, the self-heating involving the intrinsic 
graphene phonons is predicted here to lead to negative 
differential conductance. This effect has served as an 
experimental signature of the self-heating in suspended 
carbon nanotubes^. In the presence of SPP scattering 
the negative differential conductance is less pronounced 
and the self-heating effect leads to current saturation in 
Fig. [7J The different outcome of self-heating effects lead- 
ing to either negative differential conductance or current 
saturation can well be understood by the difference in 



FIG. 8: (Color online) Carrier density dependence of (a) the 
saturated velocity and (b) saturated current in graphene on 
Si02 using the self-heating model at different thermal conduc- 
tances r in kW/(K cm 2 ) from top to bottom: black - T=300 
K or r — oo, cyan - r = 4.0, magenta - r = 2.0, red - r — 1.0, 
green - r = 0.5, blue - r = 0.25, which correspond to the 
ideal thermal resistance to Si02 with an effective thickness 
h = n/r = 0, 35, 70, 140, 280, and 560 nm correspondingly. 
The dashed black curve in (a) is calculated using Eq. for 
SPP phonon in Si02- The solid squares in (b) show fitted 
saturated currents to Eq. (0 from Ref.— . 

The effect of the thermal conductance r on velocity sat- 
uration is explored in Fig. [HJ As r is reduced from infinity 
to 0.25 kW/(K cm 2 ), the saturated velocity drops by a 
factor of four. For small values of r and large densities, 
we find that Eq. ([7J) fails to fit the results of calculations 
in the full range of fields and we have chosen to fit a range 
of fields from 0.5 to 1 V//xm, where saturation is almost 
attained. The calculated saturated current in Fig. |8Jd for 
r = 0.5 kW/(K cm 2 ) agrees remarkably well with recent 
measurements in RefJ^. The data from Ref are fitted 
to Eq. ([7]) as shown in Fig. Note that accidentally 
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Eq. (|5J) reproduces fairly well the results of self-heating 
model for r = 0.25 kW/(K cm 2 ) in Fig.|5h, 
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FIG. 9: (Color online) (a) Circles show data points from 
Refi 3 . for V g =24, 21, 18, 15, 12, 9 V from top to bottom. 
The solid curves are fits to Eq. @ with r\ = 0.072 x V g 10 12 
cm" 2 , where V g is in V. The fit results are shown in Fig. [Bb. 
(b) Modeling of the current-voltage characteristics measured 
m Reh±2. The green and red dashed curves, correspondingly, 
are BTE solutions for graphene phonon scattering and both 
graphene and SPP phonon scattering at T arn b = 300 K in- 
cluding Coulomb scattering (see text). The green and red 
solid curves are self-heating model calculations with r = 0.278 
kW/(K cm 2 )^ including graphene phonon only, and both 
SPP and graphene phonon scattering, respectively. Black 
squares are experimental results from Ref.—. 

In the experiment described in Ref.—, the thermal con- 
ductance r was directly measured and does not have to be 
a fit parameter in modeling the experimental I-V curve. 
The observed mobility of about 1000 cm 2 /Vs in RefJ^ 
can be reproduced in our calculations by assuming scat- 
tering with charged impurities^ of density rii = 4.5 x 10 12 
cm" 2 and with a smaller nj = 3.5 x 10 12 cm -2 in the 
presence of SPP scattering from the SiC>2 substrate. The 
charge carrier density was fixed by the gate voltage at 



n ~ 10 13 cm -2 in RefJ^. In Fig. [9]o we show that the 
calculated current is significantly larger at high biases in 
the presence of both intrinsic graphene and SPP phonon 
scattering at room temperatures. Most importantly, the 
current does not show the saturation that is observed in 
the experiment. On the other hand, using the experi- 
mentally measured temperatures, our self-heating model 
with r = 0.278 kW/(K cm 2 ) including the SPP scattering 
very nicely reproduces the experiment. However, if we 
assume self-heating model including only the graphene 
phonons (no SPP) we do not find full saturation, even at 
fields up to 2 V/^m, and the calculation overestimates 
the measured currents at high biases. 



V. CONCLUSION 

In conclusion, we explored the effect of current and ve- 
locity saturation on the carrier density and impurity con- 
centration including intrinsic graphene phonons and SPP 
scattering on polar substrates. The observed full cur- 
rent saturation is simulated only when we account for the 
self-heating. Without self-heating, the current densities 
arc predicted to be too high for either graphene phonon 
scattering or SPP scattering. While impurity scattering 
modifies substantially the low-field mobility, it has little 
effect on the saturated velocity. Furthermore, the sat- 
urated velocity depends very weakly on carrier density 
and the choice of substrate. These dependencies served 
as the basis for invoking SPP scattering for graphene on 
SiC>2 as the velocity saturation mechanism in Ref i n i 12 . A 
more direct experimental confirmation for the SPP role 
in velocity saturation is desirable. We predict a factor of 
four enhancement of the saturation current if self-heating 
effects are minimized. This can be achieved by using an 
appropriate choice of substrate with high thermal con- 
ductivity, scaled down insulator thickness, and keeping 
the graphenc/substrate contact thermal resistance low. 
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